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1. Introduction 

Precise low-energy experiments on parity nonconservation (PNC) in atoms provide a test 
of the standard model of elementary particle interactions. In a recent PNC experiment 
with cesium by Wood et al ^ the PNC El transition amplitude between the 6s and 
7s states has been determined with an unprecedented accuracy of 0.3 %. At this level 
of accuracy a small perturbative potential (such as the Breit interaction, QED vacuum 
polarization) may influence the result. The PNC weak interaction matrix elements are 
determined by the electron wave function inside a nucleus. Having this, as well as 
some other possible applications, in mind we present in this paper a simple analytical 
approach which allows one to calculate corrections to the electron wave function near 
the origin due to a perturbative potential. As an application we consider corrections 
to weak matrix elements produced by the QED vacuum polarization (Uehling potential 
0, Wichmann-KroU potential 0]) and a QED correction to the electron potential. 
Relativistic units c = Ti = 1, = a = 1/137.04, aao = (ao is the Bohr radius) are 
used below everywhere, if not specified otherwise.. 

2. General formalism 

Our goal is to describe a variation of an electron wave function inside the nucleus due 
to some local potential V{r) considered as a perturbation. Let us presume that the 
electron motion in an atom is described with the help of a single-electron spherically 
symmetrical Hamiltonian if, which is, generally speaking, relativistic, i.e. of the Dirac 
type. The wave function for the n-th energy level ipnjiiTc) which is the Dirac 4-spinor 
characterized by the total momentum j, orbital momentum / and projection of the total 
momentum /i, the latter index is suppressed, satisfies the Schroedinger equation 

Enjllpnjlir) = Hiprijiir) , (1) 

with appropriate boundary conditions at r = and oo. The spherically symmetrical 
potential V{r) considered as a perturbation results in a variation of the wave function 
Sipnji{f) that we are looking for. This variation can obviously be presented as 

6i/jndr) = J G,i{r, r'; E^) V{r) tfj^^iir') dV , (2) 

where Gji{r, r'; En) is the corresponding Green function for the operator H, that can be 
expanded as a series over a full set of solutions of the Schroedinger equation (||) for the 
j, I wave 

G„(r.r-;E„)^i: *'f' -^y'""' . (3) 
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This Green function is necessarily a 4 x 4 matrix in the Dirac indices that are suppressed. 
Remember that the n-th electron level is excluded from summation in (^). To keep trace 
of this fact a tilde hat is used for Gji{r,r'; En) thus distinguishing it from the Green 
function Gji{r,r' , E) defined for arbitrary, non-specified energy E. Using orthogonality 
and a completeness of a full set of the wave functions i^njiiT^) that satisfy the Schroedinger 
equation (|1|), i.e. the fact that 

J^njiir)^n'jiir)d^r = Snn' , (4) 
Y: <,(r)^„,,(r') = 5(r-r') , (5) 

n,j,l,fj. 

one finds from that Gji{r, r'; En) satisfies the following non- homogeneous equation 
(En - H) Gjiiv, r'; E„) = 5{t - r') - ^„,,(r)^+,(r') , (6) 

as well as an integral condition 

J G,i{r,r';En)^nAr')d'r = . (7) 

The last term in the right-hand side of (|]) as well as condition (0) arise due to the same, 
mentioned above reason, namely that the n-th energy level is excluded from summation 
in (I). 

Our goal is to describe the wave function variation inside the nucleus. To this 
end it is sufficient to find its variation only at one point inside, for example at the 
nuclear center which we consider as an origin. Having done that one recovers the wave 
function everywhere inside (and in close vicinity of) the nucleus simply by scaling its non- 
perturbed value to comply with the variation found at the chosen point. This statement 
follows from the fact that the potentials considered in this paper are supposed to produce 
negligible effects inside the nucleus. They contribute only due to their existence in a 
nuclear exterior region. This means that equations governing the electron wave function 
inside the nucleus remain intact by the perturbation. As a result the perturbation can 
only scale the wave function inside. We will use this fact looking for the wave function 
specifically at the origin, i.e. hunting for 6ipnjiW- 

In order to find the wave function at the origin one needs to find the Green 
function G{0, r'; En)- This task can be conveniently fulfilled using the following method. 
Consider an energy E as an arbitrary parameter assuming only that E is located in 
some vicinity of a chosen atomic energy level En- Let us call ipjii'^^y E) a solution of the 
Schroedinger equation 

{E-H),p^i{v,E) = Q - (8) 

Obviously for an arbitrary energy E this solution cannot satisfy proper boundary 
conditions both at infinity and at the origin. However, we can always consider the 
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proper condition at one of the two points. We may assume therefore that ipji{r,E) 
decreases at infinity 

^ji{r,E)^0, r^oo. (9) 

Solutions of this type are necessarily singular at the origin, if the energy parameter E 
does not coincide with some physical energy level. It is convenient, nevertheless, to 
introduce the following normalization integral 

/^+,(r)^/.,,(r,E)dV = l. (10) 

which converges well because the mentioned singularity at the origin is compensated for 
by a zero-type behaviour of the regular solution. 

The function %1){y,E) possesses several important for us properties. Firstly, at 
the point E = En the boundary condition at the origin can obviously be satisfied. 
Combining this statement with equations (||),(P) and ([ToD we conclude that for E = En 
the function ipji{r, E) coincides with the wave function for the n-th energy level 

^jl{r,En) = Ipnjli^) ■ (11) 

To reveal another useful property of ipji{r,E), let us differentiate the Schroedinger 
equation (|1|) over the energy at the point E = En 

{En-H)^^t^^^^ = -^,,{r,En) , r>0. (12) 

Keeping in mind that ipji{r,E) is irregular at the origin we have to suspect that the 
right-hand side of (|T2|) may include the delta-term oc 5(r), or its derivatives. This 
suspicion, justified below, prompts to remember that r > when ( |l2l) is taken literally. 
To proceed let us differentiate the normalization condition (|10D over E at the point 

E = En 

<,W*^^dV = 0. (13) 

Compare now a set of equations (|12D and ( p!3D with (|^) taken at r = 0, r' > and (|^. 
Their obvious identity indicates that the Green function at the origin Gji{Q, r; En) can 
be presented as 

G,KO, r; En) = ^PndO)^^^%^ ■ (14) 

We can now clarify behaviour of the right-hand side of equation (|I2|) at r = 0. From 
(|I^) and (^ we deduce that our expectation was correct, at r = the right-hand side 
of ( [T^ ) should indeed be modified to include an additional delta-term ~ 6{r) (though 
for our purposes it suffices to consider this equation only at r > 0). 
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We find that iljji{r,E) is a very convenient object. It allows one to describe 
simultaneously a set of wave functions ipnji{^) (|Tl| ) as well as the Green function 
G(0, r; En) (|1^). It follows from (3) , (p!^ that the variation of the wave function due to 
perturbation V{r) can be expressed as 

5^n,m = ^PnAO) I ^^?^^^"V (r)^„,,(r) dV . (15) 

This convenient presentation is the main result of this Section. Its applications are 
discussed below. 



3. Perturbation at small distances 



Let us consider perturbative potentials V{r) located at small separations from the 
nucleus. We will assume, however, that a potential considered gives a significant 
contribution mainly outside the atomic nucleus, i.e. assume that a region of distances 
where the potential is important satisfies 

riv<r<ao/^, (16) 

where r^r is the nuclear radius. A wave function variation in the nuclear interior can be 
described by a scaling factor, which we find below. 

Conventional presentation for the Dirac four-spinor for spherically symmetrical 
potentials reads 

1'.Ar) = - ( ] . (17) 



r 



i£/n(r)%^(n) 



Here fn{f) and gni'f') are the large and small radial components of the spinor, the indices 
jl for them are suppressed. They are normalized as 

{f„,{r) + 9l{r))dr=l, (18) 

Qjifj_{n) and Qjj^{n) = —{cr ■ n) are spherical spinors, and I + 1 = 2j. The Dirac 

equation (|I|) for a motion in the potential U{r) in this notation takes the familiar form 

fLir) + -Ur) = {m + En- U{r)) g„{r) (19) 
r 

-g'Ar) + -gn{r) = {m-E^- U{r)) /„(r) , 
r 

where k = /(/ + 1)— j(j + l) — 1/4 = ±(j + l/2). For small separations of an electron from 
the atomic nucleus (|T6|) the potential U (r) can be approximated by the pure Coulomb 
potential created by the nuclear charge Z 

U{r) ~ -— . (20) 
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Additional simplification for this region comes from the fact that the energy of a valence 
electron is low and therefore for small separations (|16|) one can safely assume that 
~ m. With these simplifications the Dirac equation (19) reads 

fLir) + ^fnir)=(2m+^]g^{r) (21) 

-g'nir) + '^9n{r) = ^ 

There are two sets of solutions for these equations. One of them is regular at the origin. 
From (21) one finds that it behaves as 

U{r) = a^r^ , (22) 

Za 

Here 7 = (k^ — (Za)^)^^^. Clearly, this regular solution can be used to describe the 
electron wave function for small distances (p^Gf): fn{r) = f+{r), gn{r) = g+{r). The 
subscript + is used to distinguish this solution from the singular one. The latter, that 
will be called /„(r), (7„(r), behaves as 

/_(r)=a_r-^ (23) 

9-{9) = b-r'^, 6_ = ^ . 

This solution is used below for the description of the Green function. The explicit form 
for both sets of solutions found from (21) reads. 



r(±27 + i 

f±{r) = a± 



{8Zam) 



±7 



2(±7 - k) 



(24) 



n (r)-n + ^ ^±27 + 1) . . , . 



Here J^ix) is the Bessel function, x = y/SZamr. Numerical coeffici ents in (24), (25) are 
chosen to satisfy (22), (23). The coefficient a+ for the regular solution that represents 
the electron wave function should be found from the normalization condition for this 
function. A proper normalization of the coefficient a_ is discussed in detail below, see 

Following the approach of Section we need to replace En in (19) by an 
arbitrary value E, assuming that the corresponding solution E) with components 
f{r,E), g{r,E) behaves regularly at infinity as specified in (^. Further, we need 
to consider derivatives over the energy df{r,E)/dE = fE{r,E) and dg{r, E)/dE = 
gE{r,E). Our task is to find these functions at small r. Since f{r,E), g{r,E) are 
singular at the origin, we have to expect that fE{r,En), gE{r,En) are singular as 
well. Bearing this in mind we can neglect the regular, and therefore small, right- 
hand side of (O) when r — ^ 0. We deduce from this that for small r the functions 
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fE{r,En), gE{r,En) satisfy the homogeneous Dirac equation (19) behaving singularly 
in the vicinity of r = 0. The notation f-{r), g-{r) introduced above specifies exactly 
this solution of the Dirac equation. It follows from (23) that for small separations the 
following asymptotic conditions hold 

/s(r,E„) = /-(r) = a-r-^ (26) 
gE{r, En) = g-{r) = b_ r"^ . 

We need to continue this line of argumentation and find corrections of the order of ~ mr 
to the right-hand sides of (26). Observe firstly that according to (22) the right-hand 
side of non-homogeneous equation (p!2| ) is small enough for short distances to produce 
no corrections of the order of mr. Therefore the main correction to (26) arises from 
the homogeneous Dirac equation (21) when the mass term in the right-hand side of this 
equation is taken into account. Technically the easiest way to recover the correction due 
to the mass term is through an expansion of the explicit solutions (24), (25) in powers of 
mr. We will present the result below, in (P^, where a similar correction for the regular 
solution (22) is also included. 

To proceed we need to rewrite (^) in terms of large and small components of 
the Dirac spinor. Relations (22) ensure that variations of both components due to 
perturbation are proportional at the origin. This means that if we define the component 
ratios at the origin as a limit (5/+(0)//+(0) = [5/+(r)//_|_(r)],,.^o and 6g-^-{0)/ g-^-{0) = 
[Sg-\-{r) / g^{r)]r^o (the limit is necessary since one of the components may turn zero at 
the origin) then 

5/+(0)//+(0) = 6g40)/g40) . (27) 
Using this fact we derive from (|l^) 

^r?7^ = = r[fn{r)fE{r,E.^) + gn{r)gE{r,En)]V{r)dr (28) 

g+{0) Jo 

: U{r)f-{r) + 9+{r)g-{r) ] V{r) dr . (29) 





The last identity here arises because, as explained above, for small separations satisfying 
(p!6|) the wave function /„(r), (7„(r) and derivatives over energy fE^r, En), gE^r, En) can 
be replaced by the plus and minus solutions f±{r), g±{r) respectively. The magnitude 
of these latter solutions is governed by coefficients a± in asymptotic formulae (26). We 
need therefore to find the product a_|_a_ . This can be achieved using the following 
transformation. Multiply the non- homogeneous equation (p!2D by ip^jii^) and integrate 
over the full 3D space from which the interior of a sphere of radius e around the 
origin is excluded. Consider e > as a small parameter which is to be put to zero at the 
end of the calculations, e 0. The normalization condition for V'n(r) ensures that this 
procedure gives —1 in the right-hand side of ([T^). Integrating the Dirac Hamiltonian H 
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in the left-hand side of this equation by parts one observes that only the surface term 
sitting on the sphere survives 

9il>jiir,E) f . dllJnjlij^En) 

}e<r 




(30) 



Here we use representation of the spinors '>\}njiij) and dip{v,E)/dE in terms of their 
large and small components, compare (|1^, as well as the fact that infinity r = oo 
gives no contribution to the surface term since ipnjii''^) and dipjiir, E)/dE are regular 
there. For small radius r = e one can express ipnjii^) and dipjilr, E)/dE in terms of 
f±{r),g±{r), as was explained above, and use asymptotic formulae (22), (23), expressing 
thus the surface term via the coefficients a±,b± in the last identity in (30). We find 
from these transformations that a+6_ — 6+a_ = 1, or, remembering expressions of b± in 
terms of a± in (22), (23), find the product in question 

1 Za 

a+a_ = . 31 

2 7 

We possess now all ingredients necessary to derive the final result. Take equation (29). 
Substitute in its right-hand side expressions (24) and (25) for f±{r), gj-{r) that are 
supplemented by condition (|31|) on the coefficients a±. After that expand the resulting 
integrand that arises from the right-hand side of (29) in powers of mr. This expansion is 
both justified and necessary in view of the following reasons. The expansion is allowed 
because the perturbation V{r) is located in the region of small separations (|TB|). A 
typical radius where the potential is located is a parameter for this expansion. Necessity 
for this expansion is twofold. Firstly, the approach developed neglects the screening of 
the Coulomb field by an electron cloud, which is a good approximation only in the close 
vicinity of the nucleus. Secondly, the procedure described neglects the regular solution 
in the right-hand side of the non-homogeneous Dirac equation, which is justified only 
for small distances where this solution is small. 

Analytical calculations described above are straightforward. The final result reads 



/+(0) g+iO) 



V{r) (a + kr) dr . (32) 



Here a is a parameter with length dimension while A; is a dimensionless coefficient 

Za h , , 

a = , 33 

7 mc 

7(472 - 1) ^ ' 
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Relations (p2D,(33) are presented in absolute units, to make them more accessible for 
different applications. 

Simple formula (^) is one of the most important results of this paper. It solves 
the main problem formulated in this section presenting a variation of the wave function 
at the origin in very transparent terms, as a linear combination of the zeroth and first 
momenta of the perturbative potential. 

Note that for the short-range potentials relative corrections to the energy and wave 
function are quite different. Indeed, we can approximate the energy variation in this 
case using asymptotic relations (22) as 

6En = {njl\V\njl) ~ (a^ + bl) / V{ry^ dr . (35) 



For an arbitrary perturbation V{r) the integral in the right-hand side of this identity may 
deviate significantly from the integral in (|32D. Therefore for short-range perturbations 
the energy variation, generally speaking, cannot serve as estimate what happens with 
the wave function. 



An interesting comparison can be made with the nonrelativistic limit of (^) that 
reads 

6^ni{0) _ m 2 

Deriving this identity we use the fact that according to (33) the parameter a turns zero 
in the limit Za — *■ 0, while from (34) one derives k — >■ 2/(2/ + 1). There is a simple short 
cut derivation that leads to ( ^61) and can be used for verification of this result. Assume 
that in the nonrelativistic limit the electron motion is dominated by the kinetic term, 
which is true for short separations. Derive from this an approximation 

G,(r, r'; E^) ^ gJ^V - r'; E„) ^ , (3T) 

where Gf^\r, r'; En) is the Green function for the free motion in the l-th partial wave 
and the last identity takes into account the fact that the binding energy is negligible 
for short distances. Remembering also that for nonrelativistic motion the wave function 
behaves as ipi{r) oc r', r ^ one immediately derives (^) directly from (^ thus 
verifying relativistic equation ( p2D that we used above. To comply with absolute units 
used in (p2D we use the same units in nonrelativistic formulae (|36| ) ,(|37| ) . 

The nonrelativistic result (|36D shows that the parameter that governs variation of 
the wave function is m / V{r)rdr. This is almost an obvious result valid for a variety of 
quantum mechanical problems [Q. The relativistic result ( p2D shows that there exists 
another parameter ma JV{r)dr. It is suppressed compared with the nonrelativistic 
parameter only by a factor Za which is not small for heavy atoms. This suppression can 
be well compensated for if the potential considered increases at small separations which 
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makes a J V{r)dr larger than / V{r)rdr. In this case the found relativistic parameter 
becomes more important than the nonrelativistic one. A perturbation due to the 
QED vacuum polarization discussed in Section ^ presents an example important for 
applications. 

We can apply the results obtained above for a specific interesting case. Consider the 
parity-violating weak interaction of an atomic electron with the nucleus that mixes si/2 
and pi/2 states of an outer electron. The matrix element for this mixing (^1/21^^1-51/2) 
is saturated inside the nucleus. Therefore the variation of the matrix element for the 
weak interaction can be found simply by adding variations of Si/2 and pi/2 states given 
in (43) 

2m /■°°_. , / 2 



HPi/2 


\w 




{piM 


W\Si/2) 



2 , V{r) {a+-kr) dr . (38) 
h Jo V 3 / 

Deriving this result we take into account that essential parameters for si/2 and pi/2 
states are Ks = —1, Kp = 1, 7^ = 7p = (1 — (Za)^)^/^ = 7 and assumed that 
A; = A;, = 6/(7(47^-1)). 

4. Vacuum polarization 



Let us apply (|32D to a specific case when perturbation originates from polarization of the 
QED vacuum caused by the Coulomb field of the nucleus. In the lowest, second order 
of the QED perturbation theory this polarization is described by the Uehling potential 
i Wp(r) 

VVp(r) = - — — / exp(-2mrC) ¥{() dC , (39) 
6tc \ r J Ji 

c 

The Uehling potential (39) is singular at the origin 

Wp(r) = (—] (\n—~C-l), mr « 1 . (41) 

Sir \ r I \ mr bj 




Y(0 = ( 1 + 1 ^^^7^ • (40) 



Here C = 0.577 ... is the Euler constant. A In mr function in (|T) describes conventional 
scaling of the QED coupling constant that manifests itself for short distances. This 
scaling factor has an interesting consequence for the problem at hand. Being introduced 
in ( P^ it results in the In^ mr divergence of the integral ma J VYp{r)dr at small r. This 
divergence is eliminated by the finite nuclear size. As a result we find an estimate for 
the variation of the weak matrix element 



-2ma / VYp{r)dr a{Zaf \n^{mrN) ■ (42) 

J 3777 



{Pl/2\W\Si/2) J 3777 
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This shows that there exists the In^mrN enhancement in the problem, as was firstly 
discovered by Milstein and Sushkov using other methods. In our approach this result 
is linked with the relativistic parameter ma J V{r)dr (33) introduced in Section ^. For 
heavy atoms the In^ mr^r enhancement compensates for the additional suppressing factor 
Za in the relativistic parameter a (33) thus making this parameter dominant. 

In order to present more accurate results describing the infiuence of the vacuum 
polarization on the wave function let us substitute (39) in (p^. Integrating over r we 
find 



5U{0) Sg+{0) 2Za' 



/+(0) g40) 37r 



/oo / k \ 

iamE,{mrNC) + — \Y{() d( 

2 r°° \ 
+ Za— Ei(mr^C)nC)dC , (43) 
Sn Ji I 



_ Za^ I 3 k{2k - 1) , 2 
~ ~r V16 472 - 1 

where Y{() is defined in (40) and we used (33), (34) to present a,k explicitly. As is 
evident from (^) it is essential to take the finite size of the nucleus into account. 
We follow in (43) the simplest way cutting the divergent integral in (43) at r = tat- The 
symbol Ei{x) in (43) represents the known integral-exponent function 

Ei{x) = / exp(-xt)— . (44) 
Ji t 

Formula (43) solves the problem formulated above, giving a simple transparent 
presentation for variation of the atomic electron wave function due to vacuum 
polarization. Similarly we can find contribution of the QED vacuum polarization 
to the parity-violating weak interaction. Substituting (39) into (^) and making 
transformations similar to the ones used in (43) we find 

Hpi/2\W\si/2) _Za' [3 1 „ 4 



Za — ^"Ei(mr^C)nC)dc) • (45) 



{pi/2\W\si/2) 7 V 4 472- 1 
Equation ( ^5]) presents the weak interaction matrix element for an arbitrary atom in 



a transparent analytical form without fitting parameters. Numerical results are easily 
obtained by a straightforward one-dimensional integration in (^Sp. One only needs to 
specify the nuclear size that can be taken as rj\i = 1.2 cm where A is the 

atomic number, see P|. f For the most interesting case of the ^'^'^Cs formula (^) gives 
correction produced by the Uehling potential 0.47%. 

f Alternatively the right-hand side of ( p5[ ) can be calculated using an expansion in powers of mrN ^ 1 
that reads ^{ Za+^ {Zaf [ ( In - C - | )2 + 0.759 ] } + O(mrjv), where C ~ 0.577. This 

expansion brings ( ^5| ) to a form that is close, but not identical to the one derived in We will not 
pursue an origin for this discrepancy since calculations in the cited paper were fulfilled up to a constant 
that was eventually used as a fitting parameter. 



12 



Compare this result with other results obtained recently. Johnson, Bednyakov and 
Soff in Ref. calculated correction due to the Uehling potential for the parity- 
nonconservation in the 6s-7s amplitude in ^^■^Cs. It proves to be large 0.4%, which 
agrees with qualitative expectations expressed by Sushkov in previously. The result 
of 1^ includes, along with variation of the weak matrix element, variations of the dipole 
matrix element and the corresponding energy denominator that, combined together, 
describe a s — s mixing measured experimentally. Ref. |^ of Dzuba, Flambaum and 
Ginges confirmes this result and supplies more details providing separate variations for 
all three quantities mentioned above. It was found that variations of the dipole matrix 
element and the energy denominator, being not small, compensate each other almost 
completely. Thus the variation of the weak matrix element proves to be 0.4%. Numerical 
calculations in ^ were restricted by the logarithmic accuracy that was improved by using 
a constant as a fitting parameter to obtain 0.4% in line with 

We are interested in heavy atoms where the parameter Za is not small, therefore 
the lowest order polarization potential (Uehling potential) may be not sufficient. The 
higher order polarization potential (Wichmann-Kroll potential) was obtained in 0. To 
calculate the correction to the weak matrix element with the logarithmic accuracy it is 



enough to know this potential at small distances [10 



VwKir) = 0.092 ^ {Zaf f — ) mr < 1 . (46) 

61T \ ^ J 

The calculation with the logarithmic accuracy gives the following ratio of the Wichmann- 
Kroll correction to the Uehling correction for the weak matrix element (see (p^)): 

^^-0.184 , (47) 

SWu ln(l/(mr^)) 

For ^^^Cs this ratio is about -0.007. This confirms the statement of that the higher 
order corrections to the polarization potential are not important (this may probably be 
explained by high momenta of the electron-positron pair in the polarization loop). 

5. Large separations, non-relativistic case 



Let us apply formula ([15|) for the case when the perturbative potential V{r) is located 



in a region of distances r that satisfy the following conditions 

ao/Z <r <ao . (48) 

An example of an application here may be the calculation of QED corrections to the 
weak matrix element which originate from the atomic electron potential (see below). 

Two simplifications are possible here. Firstly, the motion can be described by non- 
relativistic equations, and, secondly, the semiclassical approach is valid here [ffl. We 
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can therefore assume that the Dirac spinors ipjiir, E), ?/'„j;(r) can be expressed in terms 
of the single-component nonrelativistic wave functions. The angular components of the 
nonrelativistic wave functions will be called ipi{r,E) and ipnii'f') respectively. Applying 
conventional semiclassical methods in the classically allowed region, which includes 
all distances specified in (^S]), we can write 



^i{r,E) = - ( 2^4^^ ^ sin ( / p{r')dr' ) , (49) 



r \^7r v{r) 

where p{r) and v{r) are a classical momentum and velocity, p(r) = mv{r), and u!{E) is 
a classical frequency 

TW^/-^L. (50) 

Generally speaking, the velocity v{r) depends on the energy E, making the period of 
the classical motion T{E) and the frequency uj{E) energy dependent as well. Recall, 
however, that we are interested in the behaviour of an outer electron whose binding 



energy is much lower than an atomic potential when r satisfies (48). This fact makes 
the velocity and momentum in the integrand in (|49|) almost independent on energy E 
in the vicinity of the n-th energy level. In contrast, oo{E) exhibits rapid variation with 
energy because the integral for the period T{E) in ( ^0]) is saturated at large distances, 
where velocity sharply depends on energy. Taking this into account we deduce from 
(^) that in the region of interest (^8|) the following equality holds 

dE -2u{E) dE ^'^^'^^ 

This shows that in the region (^) the derivative of the wave function dipi{r, E)/dE can 
be described by a simple scaling factor {l/2uj{E)){duj{E)/dE). This statement remains 
true for the Dirac spinor dipji{r, E)/dE as well because in the considered region (^) 
the spinor is proportional to the nonrelativistic wave function (|5l|). f Using now the 
fact that for shorter separations r < ao/Z the perturbation is assumed insignificant, 
we conclude that description of the perturbation by the scaling factor remains valid 
all the way down to the nucleus. This means that for all distances inside an atomic 
core < r < ao the derivative of the Dirac spinor over energy dipjiir, E)/dE remains 
proportional to the spinor itself ipji{r,E) with a scaling coefficient identical to the one 
in the right-hand side of (|1^). Using this result in (|15D we find 

6^IJndr) = p^i;nA^)5E^ (52) 



f In what follows we will not need an explicit form for relations expressing the spinor via the 
nonrelativistic wave function. A simple fact of their linear dependence will be sufficient. 
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Here uj'{E) is a shortcut notation for a derivative of oo{E) over E, while 6En is an energy 
variation due to the potential V{r) 

6E^ = {njl\V\njl) = J ^+,(r)r(r)V^„,Kr) dV . (53) 

We conclude from (|52D that the relative variation of the wave function at the origin 
is proportional to the variation of the energy level, being independent on any specific 
features of the potential. The coefficient uj'{E„) / {2uj{En)) in this formula is expressed in 
terms of the classical frequency for the electron motion. It is very simple for calculations, 
but can be simplified even further, if one needs only an estimation. Remember again 
that large separations from the atom r > oq give large contribution to the classical 
period T{E). For these distances an atomic field can be approximated by the Coulomb 
potential — e^/r created by a singly charged atomic residue. This fact allows one to 
approximate the frequency by conventional formula of Newtonian celestial mechanics 
for the Kepler problem which for the potential — e^/r read uj{E) = (7r/e^)(2i?^/m)^/^, 
where = m — E^ is the electron binding energy. Substituting this in (^) one finds 

^^uM^l^i^nM^ (54) 

where 5Bn = m — 6En is the variation of the binding energy. Thus the behaviour of 
the wave function can be described in terms of the binding energy only, f There is, of 
course, a short-cut way to derive this result. The wave function of an outer electron at 
small r is known to depend on the binding energy according to ipnii''^) — const /it', where 
n is an effective radial quantum number defined by the binding energy En = me^ / (2n^). 
Taking variation of this relation and assuming a weak influence of the perturbation on 
const one immediately reproduces (^^. 

The numerical simulation performed | shows that for the 6s 1/2 state in cesium 
atom the Kepler approximation (|5^ ) ensures an accuracy of ~ 20% for all perturbative 
potentials V{r) = const ■ exp(— pr) with 1/ao < p < 50/ao, while an accuracy of slightly 
more sophisticated formula (|52D is even higher, of the order of 2%. 

We verifled in this Section that the approach based on (|T3p gives sensible results 
(^),(^) for the region (^8|) . Now we can use this approach to estimate the influence of 
QED radiative corrections to electron-electron interaction on the weak matrix element. 
In the non-relativistic limit the Uehling potential can be replaced by a zero-range 
potential (proportional to the 5-function). A larger correction comes from the self- 
energy operator S(r, r', E) which also reduces to (5-function in the non-relativistic limit, 

f It is amusing to observe that the coefficient in the right-hand side of originates directly from 
the Kepler law that relates cubes of periods with squares of separations. 

I An atomic potential was approximated by some local potential which reasonably reproduces the 
valence electron wave function both inside and outside the atomic core. 
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i.e. S(r, r',i?) ~ 6{r — r')V^f/(r). A semiclassical formula for 6U{r) was obtained by 
Flambaum and Zelevinsky 



SUM = In -_!!^_-V'f/(r). (55) 

^ ^ 37rm2 \Up{r) - E\ ^ ^ ^ ^ 

Here U{r) is the atomic potential, while Up{r) is the atomic potential with a correction 
that takes into account the centrifugal potential which influence the p-wave electron in 
the intermediate state (see details in |11|]). As usual, this semiclassical expression is 
not valid near the turning points where Up{r) = E. However, a very weak logarithmic 
singularity does not produce any practical limitations on the applicability of (|55|). For 
the electrostatic potential V'^U{r) = —inp^r) where p(r) is the electric charge density. 
The main contribution to the Lamb shift of the energy level produced by 6U{r) is given 
by the nuclear charge (this contribution was calculated in |^). The contribution of 
the electron charge density can be calculated using the Thomas-Fermi approximation. 
The main contribution here comes from the interval ao/Z < r < ao/Z^^^. A simple 
estimate shows that the squared semiclassical electron wave function at r ~ uq/Z^^^ is 
Z times smaller than near the origin. Therefore, the electron-electron contribution to the 
Lamb shift of s-wave electron is Z times smaller than the electron-nucleus contribution. 
According to 0] the Lamb shift of the s-levels in ^^^Cs is ~ 0.1% . Then using equation 
(0) we obtain that the correction to the weak matrix element produced by the electron 
density contribution to the SU{r) is ~ 0.001%. 

6. Summary and conclusions 

Equation (^) provides a convenient framework to calculate a variation of the wave 
function of an atomic electron inside the atomic nucleus that arises due to a perturbative 
potential in an atom. Applied to the region of large distances {ao/Z < r < ao) it results 
in equation (|5^ ) and its simplified version (p4|). There is a reason that makes them 
interesting for applications for the atomic parity nonconservation. The most difficult 
and cumbersome part of theoretical investigation in the latter problem present many- 



electron correlations [T^, |T3|, |14|, |1^, j^]. The correlations take place exactly in the 
region of large distances discussed in Section ||. We can deduce from (|5^),(^) that an 
accuracy of calculations of the atomic spectrum provides a direct test for an accuracy of 
the weak matrix element calculation. There are, of course, other problems which require 
knowledge of an electron wave function in the vicinity of a nucleus such as hyperfine 
interaction, field isotopic shift, and time invariance violation. 

In the region of small separations r <C clq/Z our main result ( p2D gives simple 
transparent presentation for the infiuence of perturbation on the wave function. We 
deduce from it that for potentials singular at the origin a dimensionless parameter 
ma J V{r)dr (33) measures the strength of the potential. 
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Formula (|3^ ) was applied to find a variation of the weak electron-nucleus matrix 
element due to the vacuum polarization (Uehling potential). The result is 0.47% 
for ^^^Cs atom, which agrees with the results reported recently in P, 0- The 
contributions of the Wichmann-Kroll potential and QED corrections to the electron- 
electron interaction were found to be very small. 
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